Superfluid anisotropy in YBCO: Evidence for pair tunneling superconductivity 
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Proximity effect and pair tunneling models are applied as alternative scenarios to explain the 
recently measured ab-plane and c-axis components of the superfluid tensor in Copper-Oxide super- 
conductors which contain chains, such as YBa2Cu30e.gs- It is argued that conventional proximity 
effect models, which couple chains and planes via single electron tunneling only, are incompatible 
with the experimental observations. On the other hand several surprising features of the experi- 
mental data are readily explained by the presence of a microscopic pair tunneling process. 



PACS number: 74.50. +r, 74.20.Mn 

Introduction Recent measurements of the electromag- 
netic response of clean untwinned single crystals of 
Y Ba 2 Cu^O-;^s an d Y BaiCu^O^ reveal a large a6-plane 
anisotropy The ratio of the superfluid density in b 

and a directions (CuO chains are along the &-axis) at zero 
temperature is about rs 2.4 for the one-chain compound 
Y BCL2CU3OQ.95 and rather larger, 6 for the two-chain 
compound YBaiCu^O^ Q. Anisotropies of similar mag- 
nitudes are observed in the normal state resistivity [|3|. 
Early work on Knight shift and the NMR relaxation time 
T\ ||| in YBCI2CU3O6.95, indicates that there is an ap- 
preciable gap on the chains below T c . Taken together, 
these experiments unequivocally show that normal elec- 
trons on the CuO chains are coupled into the supercon- 
ducting state below T c . However, an unexpected feature 
of the recent data is that the temperature dependences 
of the superfluid densities in a and b directions are sim- 
ilar; both are clearly linear at low temperature |[l], and 
roughly obey p { a s) (T) / p { a s) (0) ~ p[ s] \T) / >« '(0) up to T c , 
where P^b) ^ s * ne su P ernu id density in the a(b) direction. 



Along the c-axis, pc(T) approaches its zero temperature 
value with a power higher than linear. The implications 
of these facts for the microscopic nature of c-axis coupling 
in Copper-Oxides is the subject of this paper. 

Two distinct models of plane-chain coupling are con- 
sidered here; firstly a proximity model where intrinsi- 
cally superconducting plane and chain layers are cou- 
pled through single electron tunneling only J5| and sec- 
ondly, an interlayer pair tunneling model where chains 
and planes are coupled through a Josephson-like pair tun- 
neling process. In both cases the superfluid response 
is assumed to be dominated by coherent Fermi exci- 
tations about quasi-2D chain and plane Fermi surface 
sheets ||. We mimic clean Y BCO-type electronic struc- 
tures by stacking planar "CuCV' and chain "CuO" lay- 
ers alternately along the c-axis, and employ a simple 
tight-binding description. We show that the penetration 
depth results are incompatible with the proximity model, 
but find a natural explanation within the pair tunneling 
model. 

Proximity model: The proximity model is defined by the 
mean-field Hamiltonian 



H = H + H 1 



(1) 



H v = ^2\Yl £ ^ nka c nka + e ±k (c[ k<j C2ka + h.c.) \ (2) 

k { n J 
Hl =J2 A nlnk(cl k1 cl_ kl + h.C.) (3) 
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where c\ kcr and c. 



2ka are t ne creation operators of elec- 
trons in the plane and chain bands respectively, £\k — 
— 2£(cos k a + cos kb) — p, and e 2 k = — 2i c cos fc& — p + e c . 
Here p the chemical potential and e c the relative energy 
between the chains and the planes. It is sufficient to 
assume the simplest nearest neighbor tight binding dis- 
persions for electrons on both chain and plane layers; the 
qualitative results are independent of the detailed forms 
of E\ k and 62k- ^nlnk ar e the gap functions with the 
pairing symmetries j n k on the plane {n = I) and chain 
(n = 2) bands. They are determined by self-consistent 
gap equations X n /V J2k lnk{c n k]C n ^ k {) = A„, where the 
A n 2) are the strengths of the pairing potential on 

the plane and chain layers. When A2 = 0, the chain layer 
is not intrinsically superconducting, A2 = 0, and the su- 
perconducting correlation in the chain layer is generated 
purely by the proximity effect. The coupling between the 
chains and planes in Eq. ^ is via single electron interlayer 
hopping with matrix element e± k = —2t±cos(k c /2). 
This coupling vanishes at the zone boundary k c = n. 
In the physically relevant regime the interlayer hopping 
constant t± is much smaller than the plane and chain 
hopping constants t and t c . 

The mean-field Hamiltonian (|l|) has four branches of 
quasiparticle excitations. Their energy dispersions are 



given by ±E± )k with E± yk = Jb k ± \Jb\ - c k , where 

En( A n7*fc +4fe)A and °k = (ElfcA 2 7 2fe 



-_Lfe 



b k 

£2k^ilik)' 2 + (e'ik - £ ik£2k + A 1 A 2 7ife7 2 fc) 2 . While the 
J3+,k quasiparticle excitation branch has a gap for all 
k, zero energy excitations may be present in the E— t k 
branch. The locations of such energy gap nodes are 
given by the condition Ck = 0, i.e. by the solutions of 
the simultaneous equations £ifcA27 2 fc + £2fcAi7ifc = 
and e 2 ±k - ei k e 2 k + AiA 2 7ifc7 2 fe = 0. Such solutions de- 
pend strongly on the pairing symmetries of electrons on 
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both layers. If the planar electrons have d-wave pairing 
symmetry, i.e. 71/. = cosfc& — cosfc a , line gap nodes al- 
ways exist. However, if both 71/. and 72k have isotropic 
s-wave symmetry (e.g. jik — 72fc=l) nodal lines exist 
only when the pairing functions on the plane and chain 
layers have opposite signs, A1A2 < 0, and, in addition, 
At\ > |Ai||A2|. In the somewhat artificial case where 
the pair interaction on the chain layer vanishes, A2 = 
and the "proximity" nodal line on the chain Fermi sur- 
face sheet (e2k — 0) lies along k c = tt, independent of 
the pairing symmetry of the plane layer. More gener- 
ally, proximity induced nodal line loops are present on 
the chain Fermi surface sheet. 




T/T c 



FIG. 1. Normalized superfluid densities p (s) (T)/p (s) (0) as 
functions of T/T c along a, b, and c directions in the proxim- 
ity model. 7ife = cosfc a — cosfcf, and 72ft = cosfct. t = 1, 
t c = 0.85, t± = 0.17, e c = 0.5, Ai = A 2 = 1.5 (T c ~ 0.23), and 
the electron concentration on the CuO planes is about 0.85. 

The superfluid tensor is related to the kinetic energy 
and the current-current correlation function For the 
proximity model Eq. [j], we find that the modification of 
the superfluid density along the non-chain a-axis (pa^) 
due to chain-plane hybridization is very small, and its 
temperature dependence is almost the same as for a pure 
2D system with the chains and planes decoupled. If 

(s) . 

the planar band has d-wave pairing symmetry, p a will 

be linear in T at low temperature. Along the 6-axis, 

however, the contribution of the chain band is impor- 
ts) 

tant. Close to the transition temperature, p b is dom- 
inated by the contribution of the plane band and the 
difference between pa and p^ is small. At low tem- 
perature, the leading temperature dependence of p^ is 
determined by that of the low-energy density of states 
p{u) = l/VJ2 k {S(u ~ + S(w - In the ab- 

sence of nodes, for example when 71^. = 72/t = 1 and 

(s) 

A1A2 > 0, p\ approaches to its zero-temperature value 
exponentially as the temperature decreases. Since this is 
clearly inconsistent with experiments, we consider only 
(nodeful) gapless cases. When u> < min(_E +! fc), only E-^ 



quasiparticles have contributions to p(u>). Since £7-,fc = 
y/ck/E + .k and E + ,k is weakly k-dependent around the 
gap nodes in general, p(ui) at low energy is therefore de- 
termined mainly by the structure of Ck- If A2 = 0, the 
quasiparticles around the nodal line k c = it and S2k = 
dominate the low energy excitations, and it can be shown 
that p(u>) ~ y/ui for small u. This is a consequence of the 
extreme flatness of the quasi-particle dispersion in a di- 
rection normal to the nodal line E-± ~ (k c — near 
k c = tt. Correspondingly, the low temperature superfluid 
density along the 6-direction has pronounced upward cur- 
vature, increasing rapidly as yT as zero temperature 
is approached This behavior is clearly inconsistent 
with the linear temperature dependence of the superfluid 
density in both planar directions in Y B0L2CU3O7-S and 
Y Ba,2CuiO%. On the other hand, in the realistic situa- 
tion where both Ai and A2 are finite, it can be shown 
that p(u>) ~ lo whenw« (|Ai|, |A 2 |), and hence p[ s) ~T 
when T -C (|Ai|, IA2I). In an intermediate temperature 
regime IA2I <C T <^ |Ai|, the strong upward curvature 

(s) 

in p b persists. Moreover, in the regime of physically 
relevant parameters self-consistent solutions of the gap 
equations give |A 2 | <§C |Ai|. Along the c-axis, p^ shows 
also a positive curvature. But the low temperature be- 
havior of pi s ' is more peculiar. If A 2 = 0, it can be shown 
that at low temperature p^ ~ T rather than \[T as a 
consequence of vanishing coherence factors on proximity 
nodal lines. If, however, A2 is finite and T -C A2, we find 
that p^ ~ T 2 . Figure [l] shows the typical temperature 
dependence of the three components of the superfluid 
densities in the proximity model on a normalized plot. 




Temperature 



FIG. 2. Schematic representation of the superfluid density 
in the proximity model. p^ ain and p$ ane are the super- 
fluid densities of the chains and planes when the interlayer 
coupling is zero. p$ = p^ ne + p { ^ aln . When chains and 
planes are weakly coupled by single electron tunneling, the 
chain- direct ion superfluid response develops a positive curva- 
ture near Tf ain . 
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The presence of a positive curvature in p b is not lim- 
ited to the "pure" proximity effect limit discussed above 
but is in fact a general feature in weakly coupled two-gap 
systems. Consider the case where the chain and plane 
bands are completely decoupled with transition temper- 
atures T c chain < TP lanc and superfluid densities p^iin and 
PpiLo respectively. In this case, p[ s) = p^ nc + p^L will 
have a sudden change at J^iiam ag s j 10wn m Fig ure ||. 
Switching on a weak interlayer coupling tj_ , leads to sin- 
gle transition temperature ~ TP lano but leaves a smooth 

(s) h • 

upturn in p b in the vicinity of T° 111 . 

While a d x 2_ y 2 paired state can account for the linear 

(s) 

temperature behavior of p a , we conclude that the prox- 
imity model cannot give a satisfactory explanation of the 
observed temperature dependences of p b or pi 8 '. How- 
ever, the above analysis does highlight paradoxical fea- 
tures of the experimental results, namely, the CuO chain 
layers must be intrinsically superconducting, yet there 
must be a node of the energy gap on the chain Fermi sur- 
face sheet. Otherwise, (dp ( b s) / dT)\ T = - (dpi s) /dT)| T =o 
as the chain band gives no contribution to the linear tem- 

( s) 

perature term of p b at low temperature, in contradiction 
to experiment. As a quasi- ID system has an energy gap 
under ordinary circumstances, the presence of nodes sug- 
gests that the gap function of the chain band A 2 72fc may 
have 2D character; i.e. j2k is not simply a function of fc& 
only. We are thus led to the surprising conclusion that 
superconducting coherence extends across CuO chains in 
the chain layer. 

Interlayer pair tunneling model: The failure of the single 
particle tunneling model leads us to introduce a pair tun- 
neling model of the type originally proposed by Wheatley, 
Hsu, and Anderson |J. The model consists of a kinetic 
energy term Hq, which is the same as defined in Eq. [|, 
and a local singlet pair tunneling potential term Hj: 

(4) 

where A n , r ,s = c nr fC nr+ si - c„ r xc„ r+<5T is the singlet pair 
operator. This model is compatible with the two exper- 
imental requirements specified above. (1) The pairing 
functions in the chain and plane layers must have the 
same symmetry; i.e. both s or both d. Otherwise, the 
only self-consistent solution is the trivial one A„ = 0. 
This model accounts for the presence of interchain pair- 
ing discussed above. (2) The magnitudes of the gap pa- 
rameters on chains and planes are tied together; plane 
and chain layers are not independently superconducting. 

Taking the BCS mean-field approximation, each term 
in (Q) can be decoupled, for example for the first 
term in (§), as A\ rS A^ r+ e s ~ (A^) A 2 r+f s + 

(4»,r+f *>Ai,r,« - (A\ r S )(A 2/r+hs ), where (A) is the 
thermal average of the operator A. As the system is 



translationally invariant, (A^ r s ) should be r indepen- 
dent. If we further assume the amplitude of (A' n g ) is S 
independent, then the gap function X(Ar s ) can be writ- 
ten as A(A| lr( 5) = A n a n s, here A„ and a n s are respec- 
tively the amplitude and the phase factor of X(Al irS ). 
a n s are determined by the pairing symmetries of the 
chain and plane bands. With ai$ = a.28 — o<-8, the mean- 
field decoupled pairing potential Hi can be written as 

H'j = ^\{A 2lk c\ k] c\_ ki + Ai7fc4fc T 4-*4 + h ' c ') 

k 

+^A!A 2 }, (5) 

where — J2s a $ cos ^5- On a square lattice, the 
nearest-neighbor singlet paired state can have either ex- 
tended s-wave symmetry, 7^ = cos k a + cos fcf, , or d-wave 
symmetry, 7^ = cos k a — cos fcf,. As only the d-wave pair- 
ing state can have gap nodes on the planar Fermi surface, 
we consider only this case. 

At the mean field level the proximity and pair tunnel- 
ing models H'j and Hi have the same form; the difference 
lies only in how the self-consistent procedure for the gap 
function has been applied. Thus the mean-field quasipar- 
ticle energy spectra are identical except that A171/. and 
A272/C in the proximity model spectrum should now be 
replaced by A27/C and Ait*, respectively. 

Figure || shows the components of the normalized su- 
perfluid tensor as functions of T/T c computed in the pair 
tunneling model using a typical parameter set. The ra- 
tio p b / pa at zero temperature is mainly determined 

by t c /t, while p^/pc"^ is mainly determined by t±/t. 
For the case shown in Figure H, the parameters have 
been chosen so that the ratios p b / Pa an d pi /pi at 
zero temperature are approximately 2.4 and 100, which 
gives a qualitative fit to the experimental values for 
Y Ba2CusOe.95- The overall temperature dependences 
of p[ s) (T)/p[ s) (0) and p ( a s) (T) / p { a s) (0) are quite similar. 
At low temperature p a s '(T) ~ T because of the pla- 
nar d-wave state. As both Ai and A 2 are finite in 
this case, from the previous discussion p(ui) ~ uj at 

( s) 

low energy and p b ~ T at low temperature. The 
temperature dependence of pi s '(T)/pc(0) is quite dif- 
ferent from pi s \T) I pi s \o) , especially at low tempera- 
ture where the numerical results indicates pi ~ T 2 . 
All of these results agree very well with the experimen- 
tal measurements 0j. Increasing e c so that the effec- 
tive electron concentration in the chain band decreases, 
we find that the difference between pa(T)/p a s '(0) and 
p b (T)/p b (0) becomes smaller, but the difference be- 
tween pi s) (T) / p { c s) (Q) and p { a s) (T) / p ( a s) (0) becomes even 
larger. These changes of p { a s) (T)/ p { a s) {0), p ( b s) (T) / p[ s) (0) , 
and pi^ (T) / p\f* (0) with the electron concentration in the 
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chain band agree qualitatively with the experimental re- 
sults for YBCL2CU3OQS and YBa2Cu3Oe.gr-, |J. 
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FIG. 3. Normalized superfluid densities p (s) (T) / p (s) (0) as 
functions of T/T c along a, b, and c directions in the interlayer 
pair tunneling model with d-wave pairing symmetry, t = 1, 
t e = 0.85, t± = 0.17, e c = 0.5, A = 1.5 (T c ~ 0.17), and the 
electron concentration on the CuO planes is about 0.85. 

In the above discussion we have assumed that the am- 
plitude of the order parameter (A^j) is the same in a 
and b directions. A more complete self-consistent solu- 
tion of the gap equations should allow admixture of an 
s-wave component to form a d+s wave state |10 . The 



gap nodes survive for weak admixture, but the positions 
of the nodes shift away from the diagonals of the Bril- 
louin zone. Admixture of a small s-component is con- 
sistent with the c-axis Josephson tunneling experiments 
on Y Ba2CusOT^s 0|- While full self-consistent treat- 
ment for the gap equations including a s-component is 
complicated technically in our model, preliminary calcu- 
lations suggest that this improves the agreement between 
the pair tunneling model and experiments. Another re- 
finement of the model would include a finite plane layer 
pairing interaction. 

In conclusion, we have pointed out that the planar 
anisotropy of electromagnetic response of clean YBCO- 
type structures forms a useful probe of the microscopic 
state of Copper-Oxide superconductors. We have argued 
that the data of Hardy and co-workers jij imply that the 
CuO chains layers are intrinsically superconducting but 
that there is a node on the chain Fermi surface sheet. 
The experimentally observed temperature dependences 
of superfluid densities are incompatible with conventional 
d-wave proximity models in the clean limit without an 
implausible fine tuning of parameters. The data are how- 
ever readily compatible with a model of d-wave pair tun- 
neling superconductivity. 

Acknowledgement: We thank J. Waldram for bringing our 
attention to the significance of Ref . |i| . 



[1] K. Zhang, D. A. Bonn, S. Kamal, R. Liang, D. J. Baar, 

W. N. Hardy, D. Basov, and T. Timusk, Phys. Rev. Lett. 

73, 2484 (1994); D. N. Basov, R. Liang, D. A. Bonn, W. 

N. Hardy, B. Dabrowski, M. Quijada, D. B. Tanner, J. 

P. Rice, D. M. Ginsberg, and T. Timusk, idib 74, 598 

(1995); D. A. Boon, S. Kamal, K. Zhang, R. Liang, and 

W. N. Hardy, preprint. 
[2] J. L. Tallon, C. Bernhard, U. Binninger, A. Hofer, G. 

V. M. Williams, E. J. Ansaldo, J. I. Budnick, and Ch. 

Niedermayer, Phys. Rev. Lett. 74, 1008 (1995). 
[3] T. A. Friedmann, M. W. Rabin, J. Giapintzakis, J. P. 

Rice, and D. M. Ginsberg, Phys. Rev. B 42 6217 (1990). 
[4] M. Takigawa, P. C. Hammel, R. H. Heffner, and Z. Fisk, 

Phys. Rev. B 39, 7371 (89); T. Imai, T. Shimizu, H. 

Yasuoka, Y. Ueda, and K. Kosuge, J. Phys. Soc. Jpn, 

57, 2280 (1988). 
[5] R. A. Klemm and S. H. Liu, Phys. Rev. Lett. 74, 2343 

(1995). 

[6] R. Liu, B. W. Veal, A. P. Paulikas, J. W. Downey, P. J. 
Kostic, S. Fleshier, U. Welp, C. G. Olson, X. Wu, A. J. 
Arok, and J. J. Joyce, Phys. Rev. B 46 11056 (1992); M. 
Schabel et al, unpublished. 

[7] D. J. Scalapino, S. R. White, and S. C. Zhang, Phys. 
Rev. Lett. 68, 2830 (1992). 

[8] W. A. Atkinson and J. P. Carbotte, preprint. An up- 
turn feature in p' s ' was found in numerical calculations 
of normal-superconducting layered systems by these au- 
thors. Their results are compatible with the vT depen- 
dence demonstrated here. 

[9] J. M. Wheatley, T. C. Hsu, P. W. Anderson, Phys. Rev. 
B 37, 5897 (1988); S. Chakravarty, A. Sudbo, P. W. An- 
derson, and S. Strong, Science 261, 337 (1993). 
[10] Q. P. Li, B. E. C. Koltenbah, and R. Joynt, Phys. Rev. 

B 48, 437 (1993). 
[11] A. G. Sun, D. A. Gajewski, M. B. Maple, and R. C. 
Dynes, Phys. Rev. Lett. 72, 2267 (1994). 



4 



